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1 . Introduction 

The  recent  developments  in  nonlinear 
wave  theory  over  the  last  ten  or  fifteen 
years  are  already  too  extensive  for  a de- 
tailed review  in  one  lecture.  But  it  is 
hoped,  through  a discussion  of  some  current 
ideas  and  iheir  relation  to  earlier  work,  to 
touch  upon  enough  topics  to  give  some  gen- 
eral impression  of  the  field. 

One  of  the  remarkable  advances  in  the 
area  was  the  discovery  of  highly  nontrivial 
exact  solutions  to  some  of  the  prototype 
equations.  These  were  first  found  for  the 
Korteweg-de  Vries  equation 
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was  the  solitary  wave 

r)  = So^sech^  ^ ‘ o^t).  (3) 


However,  the  unexpected  boost  came  from 
the  discovery  by  Gardner,  Greene,  Kruskal 
and  kiiura  (1967)  that  explicit  formulas 
could  be  given  for  the  interaction  of  any 
number  of  solitary  waves  (with  different 
parameters  o)  and  that  even  more  general 
solutions  could  be  related,  via  a scattering 
problem,  to  the  solution  of  a linear  integral 
equation. 

The  equation  is  similar  in  form  to 
Burgers'  equation 

(4) 


^4.  + till  = 0. 
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(1) 


This  equation  was  derived  in  1896  as  a long 
wave  approximation  for  water  waves  but  is 
of  general  interest  since  it  combines  a sim- 
ple dispersive  linear  part,  whose  dispersion 
relation  is 


w f = 0,  (2) 

with  a typical  nonlinear  term.  As  a conse- 
quence it  appears  as  an  approximation  in 
other  physical  settings:  plasma  physics  is 
an  example.  Korteweg  and  de  Vries  had 
already  noted  that  solutions  of  the  form 
17  = n(x  - Ut),  representing  wavetrains 


which  can  be  reduced  to  the  heat  equation 
(and  hence  solved  in  detail)  by  the  Cole- 
Hopf  transformation 

2^ 

u = - 2(log  ..)  = - — (5) 

x if 

But  the  methods  developed  required  a whole 
sequence  of  quite  ingenious  steps,  not  just 
a single  trick. 

perhaps  the  easiest  way  to  show  the 
surprising  novelty  of  this  work  is  to  quote 
immediately  the  result  for  N interacting 
solitary  waves.  If  the  parameters  are 
OTj Q^,  the  solution  is  given  by 


translating  with  constant  speed  U,  could  ^ 

be  found  in  terms  of  elliptic  functions  where  |D|  is 

('cnuidal  waves'),  and  that  one  limiting  case  ments 
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log|D|.  (6) 

the  determinant  with  ele- 
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Imagine  calculating  second  derivatives  of  de- 
terminants, and  ihen  even  three  more  deriv- 
atives for  n , if  a direct  verification 

XXX 

were  necessary! 

The  result  of  the  interaction  is  also 
surprising.  It  can  be  shown  from  (6)  that 
the  original  solitary  waves  eventually  emerge 
unchanged,  and  the  only  memory  of  the  inter- 
action is  a constant  displacement  of  position 
from  the  path  each  one  would  have  otherwise 
followed.  Figs.  1 and  Z are  for  other  equations 
discussed  later,  but  the  behavior  is  typical. 

After  considerable  further  development, 
with  many  people  making  contributions,  sim- 
ilar results  have  now  been  found  for  other 
prototypes  that  combine  simple  dispersion 
with  typical  nonlinear  terms.  These  include 
Modified  K.dV.:  u^  t u^u  + u "0,  (8) 

t X XXX 

Sine-Gordon:  u ^ - u + sin  u = 0,  (9) 

tt  XX  ’ ' ' 

Cubic  .Schrodinger : iu^  + u^^  *■  v|u|^u  s 0, 

flO) 


Boussinesq:  u ^ - u 

^ tt  XX 


xxxx  - 0, 

(H) 


and  certain  difference  equations.  It  is  sig- 
nificant that  these  were  all  posed  as  key 
equations  of  the  subject,  with  realistic  appli- 
cations in  mind,  before  the  possibility  of 
such  general  solutions  was  realized.  They 
were  not  manufactured  to  be  solvable. 

Z . Successive  approximations  and  perturba- 
tion exparsions. 

One  of  the  quickest  approaches  to  re- 
sults such  as  (6)  starts  in  a relatively 


straightforward  way  from  a perturbation  pro- 
cedure. Since  it  has  not  appeared  before  it 
will  be  de.scrlhod  briefly.  The  original  mo- 
tivation was  to  learn  how  some  of  the  more 
standard  perturbation  procedures,  which  have 
been  so  successful  in  other  nonlinear  prob 
lems,  would  fare  on  these  particular  equa- 
tions . 


In  the  case  of  the  Korteweg-de  Vries 
equation  (1)  it  is  convenient,  hut  not  cssen 
tial,  to  introduce 


n = - 


(!<:) 


and  take  an  equivalent  form 


1/  + 0 = 6(/-' 
t ^xxx  > 
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Then  one  would  take  the  series 


<1'  = 

1 


with  or  without  an  explicit  small  amplitude 
parameter,  and  obtain  the  hierarcliy  of  suc- 
cessive equations 

a\ 


9(//. 
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The  simplest  solution  of  (14)  is 

B(x,t)  - ye  (16) 

and  this  can  be  used  as  a preliminary  trial 
case.  After  ^ ^^<3  </')  are  calculated, 

it  becomes  apparent  and  is  readily  checked 
that 


J 


J-' 


(17) 


The  original  series  is  convergent  only  for 
large  enough  positive  x;  however  it  is 


936 


easily  sumnied  to  the  form 


*■  ■ ri- 

a 

which  is  valid  everywhere.  It  gives  in  fact 
the  solitary  wave  (3).  Notice  that  in  the 
final  result  there  is  no  limitation  to  small 
amplitude . 

Of  course  finding  the  solitary  wave  is 
no  great  achievement,  but  we  see  the  nature 
of  the  problem.  For  solutions  which  tend  tfl 
zero  as  x — oo,  we  might  start  with  a 
simple  perturbation  expansion,  but  then  its 
usefulness  will  depend  on  summing  the 
series  in  some  way  to  obtain  a form  valid 
for  all  X. 

This  approach  has  been  investigated  in 
detail  by  Rosales  (1977).  A general  solution 
of  (14)  may  be  written  formally  as  a Fouriex 
integral 

B(x)  = /e*^d\(k),  n = kx  + k’t;  (19) 


we  include  the  formal  interpretation  that  the 
integral  may  contain  contributions  from  only 
discrete  points  k = ia^  and  stands  for 

-O’  X+Q-’t 

B(x)  = 

The  t-dependence  is  not  displayed  explicitly 
for  reasons  that  appear  in  a moment.  After 
(/2  And  03  have  been  calculated  it  becomes 
clear  and  may  be  verified  in  general  that 


tl'Ax) 


i(n,  + ...+n,) 


(k,  +ic2  jfrj  + kj).  . . (kj_,+  k.) 


d\(k,  )...dX.(kj) 


(21) 
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(-D^'V  .../  B(- 

X X 

dzi  . . •dz 


*1^*.2,...B(- 
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(22) 


This  formula  replaces  the  power  in  (17)  by 

a more  general  product.  If  the  similarity 
can  be  exploited,  there  is  the  possibility  of 
summing  the  series  for  0 in  analogy  with 
(18). 

First  in  the  case  of  (20)  we  may  write 
B(^-^)  = p'^(e)p(n).  (23) 


where  p is  a column  vector  with  elements 


Pn  = ' ) 


(24) 


and  p^  denotes  the  transpose  i.e.  the  row 
vector.  Then,  (22)  may  be  written  as 

0j(x)  = (-l)^’V(x)I^''(x)P(*).  (25) 

wliere  P(x)  is  tlie  matrix 

P(x)  = j p(2)p^(z)dz,  (Z6) 


with  elements 

P_„(x)  = - — ■■  p (x)p  (x).  (27) 

mn  Q T o '^m  n 

m n 

In  (25)  we  now  have  0j  expressed  as  a 
power;  this  time  as  a matrix  instead  of  a 
scalar.  We  can  evaluate  the  sum: 

on 

1 

T -1 

= pNl  + P)  P 

= - Trace{(I  + P)’*|f) 

= - -Ij  log  det|l  + P|.  (28) 

This  gives  a form  equivalent  to  (6)  for  q. 

When  B is  not  separable  in  the  form 
(23),  the  expression  in  (22)  can  be  related 
to  the  (j  - l)th  power  of  an  operator 
applied  to  B.  To  do  this  write  B(^  ^ 
as  B(x,  y)  and  introduce  an  operator  ^ 
defined  by 

6f(x,y)  = f f(x,  z)B(z,  y)dz  . (29) 
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Then,  fronrt  the  solution  cen  be  written 

= K(x,x),  (30) 


K(x.y)  = 


Formally,  we  could  write  this  as 

K(x,y)  = (I  + B)'*B(x,y),  (32) 

but  would  then  have  to  give  a meaning  to  thf 
operator  (I  + 6)  However,  we  can 
instead  obtain  this  result  in  the  form 

(I  BlK(x,y)  = B(x,y)  (33) 

directly  from  (31),  and  interpret  it  as 


of  (15)  for  the  K.dV,  equation),  since  in 
solving  for  the  j-th  iterate  as  a nmltiple 
Fourier  IntegnU  one  divides  by 

<--(k,  ) +■...<■  u(kj)  - w(k,  f...+  k^),  (37) 

where  w = w(k)  is  the  linear  dispersion 
relation.  Most  of  this  has  to  cancel  with 
the  numerator  and  add  just  a simple  factor. 


such  us 


k.  + k. 
J-1  J 


for  the  K.dV.  equation.  At  present  the  require- 
ment of  a simple  factorization  has  limited  the 
technique  to  the  equations  mentioned.  Of  course, 
it  IS  not  necessary  to  sum  the  perturbation 


^^(x,y)  + / K(x,  z)B(z,  y)dz  = B(x,y).  (34)  series  exactly,  this  would  be  too  much  to  ex- 


This  is  the  Marchenko  integral  equation 
which  was  originally  found  as  the  end  result 
of  the  'inverse  scattering  method*. 

The  approach  goes  through  with  appro- 
priate modifications  for  the  other  equations 
noted  in  (8)- (11).  It  is  also  interesting  that 
when  it  is  applied  to  Burgers'  equation  (4) 
it  rapidly  leads  to 

j=i 

where  B is  any  solution  of  the  linear  heat 
equation.  Since  (35)  is  the  logarithmic 
series  (5),  follows  with 

^ = I + B.  (3fj 

The  key  step  is  in  finding  an  effective 
way  to  sum  the  series  corresponding  to  (21) 
for  each  case.  So  far,  this  step  relies 
very  strongly  on  a simple  factored  form  for 
the  denominator  in  the  integrand.  The  de- 
nominator arises  from  the  linear  dispersive 
operator  in  the  problem  (the  left  hand  side 


pect  n more  complicated  cases.  It  is  only  net  - 
ressary  fo  extract  the  principal  part  that  leaves 
the  remainder  uniformly  small. 

3.  Another  example  of  interacting  solitary 


Looking  beyond  these  special  equations  , 
It  is  natural  to  ask  whether  the  ret-ults  are 
typical  for  other  problems.  In  parti‘“ular, 
an  intriguing  question  is  whether  so  wy 
waves,  if  *hev  exist  for  the  problen^,  always 
interact  cleanly  and  eventually  emerge  with 
their  initial  identities  intact.  For  a general- 
ized form  of  the  Kurteweg-de  Vries  equation, 
which  we  have  studied  for  other  reasons,  the 


answer  ts  yes. 


The  equation  is 


\ ^ J V(X  - t)n^f^,t)d^  0. 


It  gives  the  Korteweg-de  Vries  equstion  when 
V(x)  = 6''(x),  (40) 

but  esn  give  sny  lincsr  dispersion  relstion 
by  chousing 
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T- 


V(x)  /'^c((()e*'‘’‘d)c.  (41) 

-00 

where  c(K)  is  the  desired  phase  velocity. 

The  original  interest  (Whitham  1967)  was  to 

soften  the  r)  term  to  allow  the  various 

XXX 

breaking  phenomena  observed  in  water  wave 
A particular  case  studied  is 

V(x)  = 2®  ®(Kl  = - j 

C v r K 

since  (39)  can  then  be  reduced  to  the  differ- 
ential equation 

^ ^ 

This  has  solitary  wave  solutions  and  a single 
one  can  be  given  analytically.  But  the 
expression  is  complicated  enough  that  it 
seems  at  present  very  unlikely  that  analytic 
formulas  could  be  found  for  the  interaction 
of  a number  of  them.  Here  the  solitary 
wave  •solution  includes  a limiting  case  of 
maximum  height  for  which  the  crest  has  a 
sharp  angle  (relevant  to  the  original  purpose 
in  water  waves). 

Dr.  Bengt  Fomberg  investigated  the 
interactions  numerically  and  found  in  every 
t.ase  that  the  solitary  waves  emerged  cleanly 
with  the  typical  phase  changes  seen  in  the 
earlier  work.  This  was  true  even  for  the 
interaction  of  the  wave  of  maximum  height 
with  a smaller  one.  See  Fig.  i. 

Other  solutions  of  (43)  were  also  com- 
puted and  a selection  will  be  shown. 

4 , The  step  problem;  the  well  problem. 

The  interaction  of  solitary  waves  has 


relate  to  the  question  of  solving  integral 
equations  such  as  (34)  when  the  function  B 
is  not  separable. 

One  basic  problem  is  to  find  the  dis- 
turbance produced  by  an  initial  step  function^ 
In  gas  dynamics  this  would  produce  a shock 
wave.  In  the  dispersive  problems  considered 
here,  the  equations  are  reversible,  there  is 
no  dissipation  and  the  solution  must  have 
quite  a different  character.  In  the  setting  of 
plasma  dynamics  such  solutions  are  some- 
times referred  to  as  'collisionless  shocks'. 

For  the  K.dV.  equation  the  fomulation 
of  this  problem  leads  to  a fairly  complicated 
B function  in  (34)  and  the  solution  of  the 
integral  equation  has  not  yet  been  found. 
However,  an  earlier  approximate  approach 
gives  the  main  parts  of  the  solution.  Apart 
from  solitary  waves,  solutions  of  dispersive 
wave  problems  tend  naturally  to  be  oscillatory 
in  character,  and  then  a detailed  description 
would  obviously  be  complicated.  With  this 
in  mind,  an  approximate  treatment  for  mod- 
ulated wavetrains  was  developed  (Whitham 
1965,  1974),  which  derives  equations  for 
various  overall  quantities  such  as  local 
amplitude  a(x,t),  local  wave  number  k(x,t)« 
and  mean  disturbance  T7(x,t).  This  theory 
was  applied  to  various  questions  concerning 
the  concept  of  nonlinear  group  velocity  and 
certain  instability  problems,  but  only  half- 
hearted attempts  were  made  on  the  step 
problem  because  of  worries  about  the  valid- 
ity of  the  averaging  procedures  near  the 
front  of  the  disturbance.  However,  Gurevich 
and  Pitaevskii  (1974)  show  that  the  theory 


naturally  drawn  a lot  of  attention,  but  there 
are  many  other  important  problems.  These 
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i^ivcs  pertectlyr  reasonable  results,  at  least 
for  the  K.dV.  equation.  In  that  case,  for  4 
step  carrying  nn  increase  in  level 

r-  U X > 0, 

n = (44) 

1.  0,  X < 0, 

the  quantities  a,  k,  rj  are  all  functions  of 
x/t,  the  similarity  form  expected  for  (44), 
and  the  disturbance  is  confined  to 

- ^ < f j . (451 


(It  should  be  remembered  that  in  the  stan- 
dard form  (1),  the  description  is  relative  to 
a moving  frame  of  reference  so  that  the 
ultimate  velocities  corresponding  to  (4S)  are 
positive).  The  solution  is  given  by 


K(8)  ^ ♦ S , 


6^‘^K(s) 


wliere  s is  the  solution  of 


2 - 

~3 — 


4 s^d  - s‘)K(s)  X 

7 E:(s)  - (1  - sMK(s)  t 


(46) 


(47) 


It  is  shown  in  Figures  3 and  4.  A typical 
result  of  direct  numerical  calculation  of  the 
K.dV.  equation  (by  Dr.  Fornberg)  is  shown 
in  Fig.  6. 

In  the  case  of  a step  that  decreases  the 
level,  i.e. 


, 0,  X > 0. 

r)  = . (48' 

- 1,  X < 0, 

the  oscillations  are  negligible  and  the  approx- 
imate solution  is  just 

n - f . - 1 < ^ < 0.  (49) 


This  is  a simple  solution  of  the  equation 

tij  ♦ = 0 (50) 

when  the  dispersive  term  neglected 

A typical  result  of  direct  numerical  calcula. 
tion  is  shown  in  Fig.  6, 


results  via  the  integral  equation  (34|;  (<19) 
looks  deceptively  easy! 

Apart  from  their  direct:  interest  the 
above  solutions  are  useful  in  bluiding  up 
more  involved  cases,  if  only  qualitatively. 
One  such  case  is  the  rectangular  well 


0, 

0 < 

X, 

- 1. 

- 1 < 

X < 

0,  (51) 

0, 

X < 

- i 

According  to 

(45)  and 

(49). 

the 

disturbance 

would  be  a conibination  of  the  results  for  th-2 
two  separate  steps  until  t - 3f  at  which 

time  the  two  parts  begin  to  interact.  After- 
wards the  solution  would  be  predicted  to  be 
something  like  the  numerical  results  shown 
in  Figs.  7 and  8. 

In  gas  dynamics  and  similar  dissipative 
systems  (also  for  Burgers'  equation  (4))  the 
asymptotic  behavior  would  just  be  the  tri- 
angular wave  shown  in  Fig.  9.  Dissipation 
would  smooth  out  the  shock  at  the  rear 
slightly  ljut  it  could  often  be  approximated 
as  a dii'continuity . In  suitable  norma''  /ed 
variables  the  approximate  solution  woul<. 
satisfy 

and  be  giver  l)y 

- s(t)  < X • 0 

>1  = (53) 

0,  X < - s(t) 

where  X - • s(t)  is  the  shock  position. 

The  determination  of  s(t)  depends  on  the 
choice  of  the  correct  shock  condition.  For 
these  familiar  cases,  it  leads  to  the  con- 
dition: 

velocity  - X strength  (54) 


i.e. 


At'*’,  strength  = At 


It  would  he  interesting  to  obtain  these 


(‘iS) 


Since  the  triangular  wave  has  width  Av 
and  depth  it  is  immediately  seen 

that  this  oosrrrsptiDd.s  to  const.;mt  area 

I Y)  dx  ~ const,  (56) 

-00 

In  these  usual  cases,  n is  proportional  to 
the  density  and  this  choice  of  shock  condition 
is  derived  from  conservation  of  mass. 

In  the  dispersive  wave  case  of  Fig.  7 , 
we  might  try  to  view  the  result  as  a trian- 
gular wave  for  the  mean  level  with  the 

oscillatory  tail  superposed  to  carry  energy 
away  (in  place  of  dissipation).  Then  we 
would  have  the  beginnings  of  a very  simple 
treatment  suitable  for  more  complicated 
problems,  A key  question  would  be  the 
'shock'  condition.  Conservation  of  mass 
would  again  suggest  (54)  and  (55).  However 
it  is  still  possible  for  mass  as  well  as 
energy  to  be  transferred  under  the  oscilla- 
tory tail.  An  alternative  suggested  by  the 
step  solution  (45)  rewritten  for  arbitrary 
step  size,  is  that  the  front  of  the  oscillatory 
part  moves  with 

velocity  “ ^ strength.  (57) 

Then 

. Is  i , 

® ” T T'  ® strength  = At  ’ . (58) 

It  is  intriguing  that  this  corresponds  to 
conservation  of 

(59) 

-00 

for  this  quantity  has  appeared  in  significant 
ways  in  other  connections.  First  i>f  all,  it 
should  be  noted  that  for  the  solitary  wave 

(3) 

Z = J (60) 

-00 


independent  of  o.  Thus  it  is  the  same  for 
all  solitary  waves,  whatever  their  amplitudes. 
Secondly,  an  imp(»rlanT  prediction  derived 
from  the  scattering  problem  referred  to  in 
Section  1,  and  not  mentioned  so  far,  concerns 
the  number  of  solitary  waves  produced  by  a 
large  initial  disturbance  r)^ix)  in  the  form 
of  a single  positive  hump  (see  Whitham  1974). 
The  number  is 

N = — — f no  dx.  (61) 

It  is  interesting  that 

NZ  = J (62) 

-00 

If  the  quantity 

r®  V> 

f dx,  (63) 

-00 

which  looks  like  some  kind  of  'action',  were 
conserved,  the  resulting  solitary  waves 
would  carry  it  all.  Somehow  they  have 
picked  up  more  before  emerging.  Rut  the 
amount  is  still  proportional  to  the  initial 
input. 

The  relation  of  these  results  with  (58)-' 
(59)  is  still  obscure;  it  is  yet  one  more 
instance  of  the  cross-relations  that  hint 
there  is  much  more  than  we  know. 

Dr.  Fornberg's  numerical  calculations 
indicate  that  the  width  of  the  triangular  wave 

i 

is  proportional  to  t^  as  in  (58)i  but  the 

_z 

strength  decays  more  slowly  than  t~^, 

closer  to  t'^  in  fact.  It  is  also  interesting 

that  the  results  of  water  wave  experiments 

-i 

by  Dr.  J.  Hammack  came  out  with  t ^ for 
the  strength;  any  viscous  corrections  would 
change  the  power  even  further  away  from  Z/3. 
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The  rectangular  well  problem  has  been 
discussed  on  the  basis  of  the  integral 
equation  (^4),  A first  version  by  Ablowitz 
and  Newell  (1973)  was  incorrect:  a new 
version  by  Ablowitz  and  Segur  (1977)  appar- 
ently favors  the  power  laws  in  (58). 

5.  Wavetrain  instabilities. 

The  modulation  theory  for  oscillatory 
solutions  is  based  on  the  existence  of  exact 
solutions  that  represent  periodic  wavetrains. 
In  some  cases  these  wavetrains  are  unstable 
to  the  modulations.  This  means  that  small 
perturbations  grow,  but  there  is  no  impli- 
cation that  the  further  behavior  is  chaotic. 
For  the  simplest  problems,  the  modulation 
equations  for  the  amplitude  a(x,t)  and 
wavenumber  k(x,t)  are 

kt  + = 0,  (64) 

(a')^  + 

where  the  frequency  u>  is  given  by 

u = uo(k)  + u,j(k)a^  - (<>6) 


u»o  (k)  is  the  linear  dispersion  relation  and 
u»2(k)  is  the  coefficient  of  the  nonlinear 
correction.  For  the  uniform  wavetrain  a 
and  k are  constants,  and  it  is  easily 
shown  that  perturbed  solutions  in  the  form 

6a,  6k  ac  ((,7, 


have 


V = ± (0)0  0)2*^ 


(68) 


If  W0W2  > 0,  the  values  of  V are  real 
and  give  nonlinear  generalizations  of  the 
linear  group  velocity.  If  ujoW2  < 0,  small 
modulations  will  grow  when  /u  is  in  the 
range 


0 < ^^<  4 I a-  . 
The  growth  rate  depends  on 


(b9) 


#i(lu)oW2  ^ (70) 

and  is  maxinmm  for 

(71) 

These  results  can  also  be  derived 
froit^  a discus  ion  of  'side  band'  interactions, 
following  Renjaniin's  original  analysis  of  deep 
water  waves.  This  will  be  indicated  for  one 
of  the  simplest  examples: 

u t 3u^u  * u - 0.  (7i!) 

t X xxx 

Solutions  of  the  form 


(t)e 


ikx 


JA, 


(t)e 


i(k-ji)x 


f complex  conjugates,  (73) 

are  considered,  where  (he  'side  bai  s' 

Ai,A2  are  initially  small.  In  the  non- 
linear interactions,  various  products  feed 
back  on  to  ihe  original  components  ano  when 
these  are  IniTuckni  hav<* 

i If”  ^ + 4a,  A, 

+ -^A.AjlAo  t “AiA^Aq, 

(74) 

TThi  ' (k-n)^A,  I^AoAo  ^ 4^iA, 

* jAjAjjA,  * 

(75) 


with  . similar  equation  for  Aj  . If  these 
equations  are  linearized  on  the  (jasis  that 

‘'o* 

the  results  in  (69)-(71)  are  recovered,  in 
this  case  with 
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Uq  = - t 


u,  = |k.  (76) 

The  growth  rate  in  (70)  is 

3kn(|'  - ^‘y\  (77) 

and  IS  a maximum  when 

M = a/2.  (78) 

As  the  side  bands  grow  at  the  expense 
of  the  main  wave,  the  linearized  approxi- 
mation ceases  to  apply  and  the  full  equation 
(74)-(75)  must  be  used.  This  is  a conser- 
vative system  and  the  energy  oscillates 
between  the  three  modes.  Eventually  the 
side- bands  decay  again  and  return  the  energy 
to  the  main  wave.  In  detail,  it  can  be 
shown  that 

iAkJi  + = constant,  (79) 

k k- M k+p 

I - ^ = constant,  (80) 

k-^  ' 

and  an  equation  for  |A]  say,  can  be 

obtained  which  has  solutions  in  elliptic 

functions.  One  limiting  case  is 

k~  ijt  k+ji  2 k 

= 

(12*i^  + a^  )+(l6fi^ -a^  )cosh  6k^(^  -^^)'^t 

(81) 

This  shows  the  growth  of  the  side  bands  in 
accordance  with  (77),  but  a maximum  is 
reached  and  the  side  bands  decay  again. 

The  other  solutions  are  a periodic  repetition 
of  this  type  of  behavior. 

For  the  sidebands  of  maximum  growth 
rate,  ^ = a/2,  the  maximum  value  in  (81) 
gives 

|A„r=4a^  |A.  = 0 -£)4a^ 

|Aj  1^  = (1  + (82J 


The  time  taken  from  10%  of  maximum  to  the 
maximum  is 

T = 2.5/ka*  (83) 

Numerical  studies  of  this  problem  were 
made  by  Fornberg  by  direct  computation  on 
(72)  and  the  results  are  extremely  interesting 
The  exact  periodic  wavetrain  was  introduced 
initially,  and  for  some  time  the  disturbance 
just  translated  in  space  without  change  of 
form,  as  it  should  ideally.  However,  the 
small  numerical  errors  were  eventually 
sufficient  to  trigger  the  bursts  shown  in 
Fig.  10.  The  clean  return  to  the  periodic 
wavetrain  (which  was  repeated  as  long  as 
the  calculations  were  run)  fits  the  above 
analysis.  To  compare  in  more  detail,  the 
various  harmonics  in  the  solution  were 
analyzed  and  a typical  case  is  shown  in 
Fig.  11.  In  these  computations,  the  problem 
is  formulated  to  be  periodic  over  a large 
space  interval;  as  a consequence,  the  wave 
numbers  present  have  to  be  integer  multiples 
of  a basic  unit.  For  the  case  shown  in 
Fig.  10  and  11  this  unit  is  ff/64.  The 
initial  wavetrain  has  a = 0.2,  k = 7 units. 
According  to  (78), 

^ = 0.1^2  units, 

for  the  most  rapidly  growing  sidebands.  In 
Fig.  11,  the  strong  growth  of  side  bands  S 
and  9 agrees  qualitatively  with  this  prediction 
But  other  side  bands  k ± 2^,  k ± 3^,  etc^  ^ 
are  also  produced  by  nonlinear  interactions 
which  were  neglected  in  the  three  mode 
equations  (74)- (75).  For  this  case  the 
values  in  (82)  are 

|Ao  I = 0.13.  |A,  I = 0.09,  |Aj  | = 0.12 
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The  extra  drain  on  |Ao  | can  be  attributed 
to  loss  of  energy  to  the  other  side  bands. 


The  time  scale  (83)  is  T - 180,  which 
checks  quite  well. 

The  bursts  repeat  with  more  and  more 
sidebands  called  into  play  at  each  burst,  and 
slowly  more  activity  in  between.  The  case 
in  Figs.  10  and  11  was  run  to  t - 7600 
with  ten  bursts  recorded. 

This  research  was  supported  by  the 
Office  ot  Naval  Research,  U.S.  Navy. 
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Fig.  1.  Interacting  solitary  waves  for  eq . (8) 


Fig.  2.  Interacting  solitary  waves  for  eq . (43) 
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